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Abstract 

We report the first observation of S° Dt'^~ and an improved measurement of Dt^~ 
based on 274 x 10^ BB events collected with the Belle detector at KEKB. We measure the branching 
fractions B[B^ D+K') = (2.93 ± 0.55 ± 0.79) x lO^^ and B{B^ -D+tt") = (1.94 ± 0.47 ± 
0.52) X 10"^ 

PACS numbers: 13.25.Hw, 14.40.Nd 
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The unitarity of the Cabbibo-Kobayashi-Maskawa (CKM) matrix P| is a crucial compo- 
nent of the Standard Model, one that is currently being tested in multiple aspects through 
B meson decays at the 5-factories. The decays -D^vr" are expected to proceed dom- 

inantly through a spectator process involving the h ^ u transition and are thus considered 
a promising mode for the precise determination of \Vub\ [2]. The decay B^ D^K~ is not 
directly accessible through a spectator process, requiring either final state interactions or 
a non-spectator VT-exchange [3|. The decay B*^ — >■ D^K~ was observed by Belle and 
confirmed by BaBar j^. For B^ —> D^7i~, both groups reported evidence 0,1^. 

In this Letter we report improved measurements of B^ —>■ DfK~ and B^ -D^tt" decays 
with the Belle detector |^ at the KEKB asymmetric energy e^e~ collider f^. The results are 
based on a data sample, collected at the center-of-mass (CM) energy of the T{4S) resonance, 
which contains 274 x 10^ produced BB pairs. 

The Belle detector is a large- solid- angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere 0. Two different inner detector 
configurations were used. For the first sample of 152 million BB pairs, a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used; for the latter 122 million BB pairs, 
a 1.5 cm radius beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were 
usedH. 

Charged tracks are selected with requirements based on the average hit residual and im- 
pact parameter relative to the interaction point (IP). We also require that the transverse 
momentum of the tracks be greater than 0.1 GeV/c in order to reduce the low momentum 
combinatorial background. For charged particle identification (PID) the combined infor- 
mation from specific ionization in the central drift chamber (dE/dx), time-of-fiight scintil- 
lation counters (TOF) and aerogel Cerenkov counters (ACC) is used. At large momenta 
(> 2.5 GeV/c) only the ACC and dE/dx are used. Charged kaons are selected with PID 
criteria that have an efficiency of 88%, a pion misidentification probability of 8%, and neg- 
ligible contamination from protons. The criteria for charged pions have an efficiency of 89% 
and a kaon misidentification probability of 9%. All tracks that are positively identified as 
electrons are rejected. 

Neutral kaons are reconstructed via the decay Kg tt+tt^. The two-pion invariant mass 
is required to be within 6 MeV/c^ (~ 2.50") of the nominal mass, and the displacement 
of the TT+vr" vertex from the IP in the transverse r-0 plane is required to be between 0.1 cm 
and 20 cm. The directions in the r-0 projection of the candidate's fiight path and 
momentum are required to agree within 0.2 radians. 

We reconstruct mesons in the channels — » 07r+, K*^K~^, and KgK~^ (inclusion 
of charge conjugate states is implicit throughout this report). (K*^) mesons are formed 
from the K^K~ {K^ti~ ) pairs with invariant mass within 10 MeV/c^ (50 MeV/c^) of 
the nominal (f) {K*^) mass. We select mesons in a wide (±0.5 GeV/c^) window, for 
subsequent studies; the M{Ds) signal region is defined to be within 12 MeV/c^ (~ 2.5o-) of 
the nominal mass. candidates are combined with a charged kaon or pion to form a B 
meson. Candidate events are identified by their CM energy difference, AE = - Ei) — Eb, 
and the beam constrained mass, Mbc = ^/E'^ — {J^iPiV^ where Eb = v^/2 is the beam 
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energy and pi and Ei are the momenta and energies of the decay products of the B meson 
in the CM frame. We select events with Mbc > 5.2 GeV/c^ and |AE| < 0.2 GeV and define 
the B signal region to be 5.272 GeV/c^ < Mbc < 5.288 GeV/c^ and |AE| < 0.03 GeV. The 
Mbc sideband is defined as 5.20 GeV/c^ < Mbc < 5.26 GeV/c^. We use a Monte Carlo (MC) 
simulation to determine the efficiency j^. 

To suppress the large combinatorial background that is dominated by the two-jet-like 
e^e^ qq {q = u, d, s and c quarks) continuum process, variables that characterize the 
event topology are used. We require | cos6'thr| < 0.80, where 6'thr is the angle between the 
thrust axis of the B candidate and that of the rest of the event. This requirement eliminates 
77% of the continuum background and retains 78% of the signal events. We also define 
a Fisher discriminant, JF, which is based on the production angle of the B candidate, the 
angle of the B candidate thrust axis with respect to the beam axis, and nine parameters 
that characterize the momentum flow in the event relative to the B candidate thrust axis 
in the CM frame 0]. We impose a requirement on JF that rejects 50% of the remaining 
continuum background and retains 92% of the signal. 

We also consider possible backgrounds from qq events containing real mesons. These 
events peak in the M{Ds) spectra but not in the AE and Mbc distributions. Thus it can be 
eliminated by fitting the AE distribution. 

Other B decays, such as B^ D~^7t~, K~7t~^7i^, with one pion misidentified as 

a kaon, require particular attention because they have large branching fractions and can 
peak in the Mbc signal region. The reconstructed invariant mass spectra for these events 
overlap with the signal mass region, while their AE distribution is shifted by about 
50 MeV/c^. To suppress this background, we exclude event candidates that are consistent 
with the K~ti^tt^ mass hypothesis within 15 MeV/c^ (~ 3o") when the two same- 

sign particles are considered to be pions, independently of their PID information. For the 
K^K^ mode there is a similar background from B^ D^n , Kgii'^. In 

this case we exclude candidates consistent within 20 MeV/c^ (~ 3cr) with the —>■ K^n^ 
hypothesis. 

Possible backgrounds from B decays via b c transitions {B — > DgDX) are also con- 
sidered. The Dg mesons from these decays have a lower momentum and are kinematically 
separated from the signal. We analyzed a large MC sample of generic BB events and found 
no peaking backgrounds. 

Another potential BB background is charmless B^ — ^ K~K^K^7T^{KgK^K^). Such 
events peak in the AE and Mbc spectra, but not in the M(D,,) distributions. They tend to be 
dominated by quasi-two-body decay channels such as (f)K^'*^'^ To reduce this background, 
we reject events with low (< 2 GeV/c^) two particle invariant masses: Mx-tt+ and M^k- for 
the Df 07r+ channel, Mk+k- and Mji*OK- for K*^K~^, and Mk+k- and M^o^- 

for — >■ KgK~^. The remaining background from these sources, if any, is excluded by 
fitting the M{Ds) distribution. 

The scatter plots in AE and Mbc for the B^ D^K~ and B^ D^t^^ candidates 
in the M{Ds) signal region are shown in Fig. P a significant enhancement in the B signal 
region is observed. Figure El shows the M{Ds) spectra for selected B^ — > DfK~ and 
B^ D~^TT~ candidates in the B signal region. In addition to clear signals at the mass 
in Fig. 121 we observe peaks at the mass, corresponding to B^ — > D~^7^~ and B^ D~^K~ , 

Our studies have shown that the backgrounds may peak in the signal region of M{Ds) 
or of AE (and M^c) but not in both simultaneously. To extract our signal, we therefore 
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FIG. 1: AE versus Mbc scatter plot for the S° D+R- (top) and B° D+tt' (bottom) 
candidates in the M{Ds) signal region. The points represent the experimental data and the boxes 
show the B meson signal region. 




M(Dg) (GeV/c^) 

FIG. 2: M{Ds) spectra for S° DfR- (top) and S° DfTr' (bottom) in the B signal 
region. The points with errors represent experimental data and the curves display the results 
of the simultaneous fit described in the text. The small feature near 1.92 GeV/c^ arises due to 
favoured decays with pions misidentified as kaons. 
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TABLE I: Results on the signal yields and branching fractions. Decay channel, signal yield from 
two dimensional M{Ds) - AE, one dimensional M{Ds) and /S.E fits, detection efficiency branching 
fraction and statistical significance. The efficiencies include intermediate branching fractions. 



Mode 








M{Ds) - AE 


M{Ds) 




e, 10-3 


B (10-5) 


Significance 


B0-> 


D+K- 
D+K- 
DfK- 


,D+- 
,Dt- 


K*"K+ 


18.2^5-0 
13.11^^ 

o '7+2.8 
' -2.2 


11-21^6 
3.6l|? 


is.itl:^ 

13.9ti^9 


2.05 
1.49 
0.86 


3.25l^;?? ± 0.88 
3.2lj;J]!^±0.87 
1.57lJJ0±0.42 


6.2(7 
4.6f7 
l.Sfj 




D+K- 


, simultaneous fit 




32.9i|^ 


36.7±^i 


4.29 


2.93 ±0.55 ±0.79 


7.9cr 




D+t:-. 












2.29 


2.03lS]:^^ ± 0.55 


4.5f7 



B° ^ Dfn-, Df ^ K*°K+ 7.0tti 5-3±3.2 ^-^-tl 1-65 1.54^°;^^ _^ 0.42 2.1cr 
fiO ^ Li+TT-, D+ ^ 5.8l^;7 6.lt^;g 6.lt^j 0.89 2.38li;i^ ± 0.64 2.4cr 

5° ^ D+TT", simultaneous fit 25.71^ ;^ 26.3t^;7 27.31^;^ 4.83 1.94 ± 0.47 ± 0.52 5.5a 



perform a binned maximum likelihood fit to the two-dimensional distribution of data in 
M{Ds) and AE, separating the backgrounds from the signal component, which peaks in 
both. We select events from Mbc signal region for this fit. For each of the three Dj" decay 
channels the AE range, -0.1 GeV< AE < 0.2 GeV, is divided into 30 bins and the M{Ds) 
range, 1.5 GeV/c^ < M{Ds) < 2.5 GeV/c^, into 200 bins. All bins in all Df submodes are 
fitted simultaneously to a sum of signal and background shapes. The signal is described 
by a two-dimensional Gaussian, with widths in both dimensions obtained and fixed using 
reconstructed signals in the data from D^7i^{D^ K^ti^ti^ , Kstt^)- The signal 

amplitude is constrained to correspond to the same branching fraction B{B^ — > Dfh~) for 
all three submodes. The fit also includes an additional two-dimensional Gaussian for 
D+h- decays. 

The background includes three components: combinatorial (linear in M{Ds) and AE), 
qq events that peak in M{Ds) and are flat in AE, and B decays that peak in AE and 
are flat in M{Ds). The levels of the three components are allowed to vary independently 
in the three reconstructed Df modes. The fit results are given in Table HI The statistical 
significance quoted in Table HI is defined as ln(£o / ^max) , where Cmax and Cq denote 
the maximum likelihood with the fitted signal yield and with the signal yield fixed to zero, 
respectively. The results of one-dimensional fits to the M{Ds) and AE distributions are 
also shown in TableUfor comparison. Figures El and El show the M{Ds) and AE projections 
for events from the signal region. The sum of the fitted signal plus background is shown by 
the solid lines while the background shape including the peaking background is indicated 
by dashed lines. The peaking background is found to be 3.6 ± 0.6 and 3.4 ± 1.3 events for 
S° D+K- and B° D+ti', respectively. 

B^ D*^h~ final states, where the low energy photon from the D* D^^ decay is 
missed, can populate the B^ D^h" signal region. These would produce a long tail on 
the negative side of the AE distribution. In theoretical models based on factorization, the 
B^ — * D*"^7r~ and B*^ —>■ -D^vr" decay widths are predicted to be approximately equal; 
there are, however, no corresponding predictions for B^ — >■ Di*^^K~ decays. To study the 
sensitivity of the measured branching fraction to a possible B^ D*'^h~ contribution, we 
perform a fit with an additional B^ D*^h^ component included, where the signal shape 
is fixed from the MC and the branching fraction is left as a free parameter. The resulting 
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FIG. 3: AE spectra for DfK^ (top) and Dfn^ (bottom) in the B signal region. 

The points with errors are experimental data and the curves are the results of the simultaneous fit 
described in the text. 



2% difference in the —>■ D^tt^ event yield (compared to the results presented in Table H} 
is added to the systematic uncertainty; the change in the DfK~ yield is less than 

1%. We also check for crossfeed between B^ —>■ DfK~ and B^ —>■ Dfn^ due to kaon/pion 
misidentification. To study this we include the crossfeed contributions in the simultaneous 
fit, with shapes fixed from the MC and misidentification probabilities obtained from data; 
the uncertainty due to this effect is found to be negligible (< 1%). 

As a check, we apply the same procedure to B^ — > D~^n~ and B^ D^K^, 
07r+,i?*Oir+,irOfs:+, and obtain B(B° ^ D+n-) = (2.8±0.1) x 10"=^ and i3(i?° ^ D+R-) = 
(2.4 ± 0.4) X 10-^ [12], which agree well with the PDG values 

The following sources of systematic error are found to be the most significant: tracking 
efficiency (1-2% per track), charged hadron identification efficiency (1% per particle). Kg 
reconstruction efficiency (6%), signal-shape parameterization (5%) and MC statistics (2%). 
The tracking efficiency error is estimated using t] decays to 77 and 71^ 71^ 71^. The K/tc 
identification uncertainty is determined from D*~^ D'^tc'^, — > K^tt^ decays. We 
assume equal production of B^B~ and B^B^ pairs but do not include an additional error 
from this assumption. The uncertainty in the meson branching fractions, which is 
dominated by the 25% error on B{Df 07r+), is also taken into account. The overall 
systematic uncertainty is 27%. 

In summary, we report the first observation of B^ —>■ Dfir^ and an improved measurement 
of ^ DfK-. We find B{B'^ DfR-) = (2.93 ± 0.55 ± 0.79) x 10"^ and B{B° 
D^Ti^) = (1.94 ± 0.47 ± 0.52) x 10~^. Since the dominant systematic uncertainty in both 
measurements is due to the branching fraction for —>■ (pn^, 3^-,^, we also report B{B^ 
DfK-)xB^^ = (10.5±2.0±1.0)xl0-^ and B^B^^ ^ Dtn-)xB^^ = (7.0±1.7±0.7) x 10^^ 
These results are consistent with previous results 0, 1^ and have better accuracy. They 



8 



supersede results from j^. 
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